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Abstract 

Using the nonlinear constraint and Darboux transformation methods, 
the (mi, • • • , mjy) localized solitons of the hyperbolic su(N) AKNS system 
are constructed. Here "hyperbolic su(N)" means that the first part of 
the Lax pair is \P y = J& x + U(x, y, t)\P where J is constant real diagonal 
and U* = —U. When different solitons move in different velocities, each 
component Uij of the solution U has at most mirrij peaks as t — * oo. This 
corresponds to the (M, N) solitons for the DSI equation. When all the 
solitons move in the same velocity, U^ still has at most mirrij peaks if the 
phase differences are large enough. 

1 Introduction 

Since the discovery of localized solitons for the DSI equation, §1 the (M, N) 
solitons, il especially the (N,N) solitons. are discussed in various ways, such 
as the inverse scattering method, iiEEMl], the binary Darboux transformation 
method Jp the Hirota method etc.O On the other hand, the nonlinear constraint 
methodB has been developed for 1+2 dimensional problems since the work of 
12] , |[]. Using the Darboux transformation in higher dimensionsBil, the solutions 
can be obtained explicitly. In [[T/J, the non- localized solutions of the N-wave 
equation were derived in this way. In | 18| , the localized soliton solutions of the 



hyperbolic su(N) AKNS system were discussed by nonlinear constraint method 
which transforms the original problem to a 2N x 2N linear system which separates 
all the variables. It shows that, if the velocities of the solitons are all different, 
there are Darboux transformations of rth degree such that the asymptotic solution 
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has at most r 2 peaks. However, the velocities are not arbitrary. For example, the 
velocities of "1-soliton" should be along the y-axis. In [Of, a localized one-soliton 
solution of the DS1I equation was obtained by nonlinear constraint method which 
transforms the original problem to an (N + 1) x (iV + 1) linear system which 
separates all the variables. 

In this paper, we show that under the nonlinear constraint, there is a way to 
construct Darboux transformation so that the number of peaks and the values of 
velocities have more freedom. This construction is valid for all the lower ordered 
equations in the hyperbolic su(N) AKNS system, including the DSI equation, 
the N-wave equation etc. For the DSI equation, this corresponds to the (M, N) 
solitons. 

The main conclusions are: (1) There are m = (mi, • • • , mjv) solitons 
which are localized (tend to as (x,y) — > oo); (2) if different solitons move in 
different velocities (characterized by an algebraic condition (|47|) and hold for the 

DSI equation), then the component of has at most mjm k peaks as t — ► oo; 

(3) If all the solitons move in the same velocity (hold for the N-wave equation), 

then zffr also has at most mjm k peaks as the phase differences are large enough. 



2 The system and its nonlinear constraints 

Here we consider the hyperbolic su(N) AKNS system 

& v = M x + U(x,y,t)& 

3=0 

where d = d/dx, J = diag(Ji, ■ • • , Jjy) is a real constant diagonal N x N matrix 
with mutually different entries. U(x,y,t) is off-diagonal with U* = —U. In this 
case, we call (|I|) a hyperbolic su(N) AKNS system, since J is real and U G su(N). 
The integrability conditions of (|I|) are 

vAi\ = K - JV A - % n A + E c::i(v k &i- k u) A (2) 

k=0 

V° - JV 3 P X = [U, V 3 A ] D - £ C:ii(V k V- k U) D (3) 

k=0 

n-1 

Ut = V n A y - JV n A x - [U, V n ] A + Y: (V k d n ~ k U) A (4) 

k=0 

where the superscripts D and A refer to the diagonal and off-diagonal parts of a 
matrix. As mentioned in jill , (|3|) and (|J) give a system of nonlinear PDEs of U 
and all the diagonal parts of Vj s, where V^-'s are determined by (0) respectively. 
Usually, only U is physically significant. 



2 



There is a connection of (|1|) with the linear system 

fiXI iF\ fiXJ+U UF 

x ~ \iF* ) ' v ~ V iF*J 

where F, Wj, X j} Zj are NxK,NxN,NxK,KxK matrices respectively 
(K > 1) and satisfy W* = —Wj, Z* = —Zj. This is a slight generalization of the 
linear system in |18j, where K should be N. 



The integrability conditions of (||) include 

(6) 



F y = JF X + UF 



iF t = X n>x + %W n F - iFZ n 

iX j+1 = X jiX + iWjF - iFZj 
W j>x = -iFX* - iXjF* 

Z jtX = iF*Xj + iX]F (7) 
z[J, W j+ i] = W hy - [U, W 3 ] + iJFX* + iXjF*J 
Z j:V = il'JXj + iX*JF 

U X = [J,FF*} (8) 
U t = W n , y - [U, W n ] + iJFX* n + iX n F*J (9) 

For U = 0, F = 0, (0) implies that Wj(X) = iQ^t), Xj = 0, Zj = iZ](t) where 
i?j(t)'s are real diagonal matrices and Z°(t)'s are real matrices. 

When Zj(t) = Q(t)lK {Ik is the K x K identity matrix) where (j{t) is a real 
function of t, (@) is just the Lax pair (JJ) for n = 1,2,3. (0) and (|) give the 
recursion relations to determine Wj, Xj, Zj, together with the evolution equations 
corresponding to (Q)-(ffl), which are the integrability conditions of (ffj). (^|) gives a 
nonlinear constraint between U and F. 

This system includes the DSI equation and the N-wave equation as special 
cases. 

In order to consider the asymptotic behavior of the solution U, here we 
suppose Qj is independent of t and Q = 0. Denote Q = J2]=o J2jA n ~ J and write 

Q = diag(a>i, ■ • • , ujn)- 

We need the following symbols and simple conclusions. 

If Mi, M.2 are j x k matrices, we write Mi = M2 if there is a nondegenerate 
diagonal k x k matrix A such that M 2 = Mi A. 

If L is a diagonal matrix, then Mi = M2 and detMi 7^ imply M\ I.M { ; = 
M 2 LM 2 -\ 

Let 

»-{: T) ^ 
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^=( A ° \ r..)- (ID 



where v 7^ is an K x 1 vector, a 7^ is a number. Let 
Then we have 

A \ -av AI K - vv* J v ' 

MAM -i = W^ + (Ao-A )|a| 2 (Ao-A„)a^ \ 

V (A - X Q )av X AI K + (A - X )vv* J v ; 

where A = v*v + |a| 2 . Moreover, 

lim MAM- 1 = ( Ao T r ^ (14) 

, ^° \ 

lim M AM' 1 = -r T .vv* (15) 

«-o I A /x + (A -A )— ) v 7 

For an (N + K) x (N + K) matrix M, denote M Bn be an N x N matrix 
containing the first N columns and N rows of M. 



3 Darboux transformation and soliton solutions 

Now we construct the solutions from U = 0, F = 0. A procedure to construct 
the Darboux transformation was proposed in . 

Let A Q (a = 1, 2, • • • , r) be r complex numbers such that A a 7^ A^ for a 7^ (3 
and X a 7^ A/3 for all a, (3. Take 

A a = diag(A a , ■ • ■ , A Q , X a , ■ • ■ , X a ) (16) 

-v ' " v ' 

N K 

Ha= l e MQa(s)) -exp(-Q a (s)*)C*\ ^ 
C a Ik J 



where 



Qa(s) 



a a j, b a j are constants, 



'a al s + b al \ 
a a NS + b aN J 



C a = (0,...,0,K a ,0,...,0) (19) 



constant x 1 nonzero vector which appears at the Z a 's column of C a . 

Let 

p a = Re(a a: iJ (p a = Im(a aIa ) vr Q = b a , (20) 
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The Darboux matrices for such {A a , H Q } can be constructed as follows. Let 

G W (X) = X-H 1 A 1 H{ 1 = G {l \\ a )H a (a = 2, 3, • • • , r) 

G^(X) = X-HPa 2 hP- 1 HV = GV\K)HV (a = 3, 4, ■••,*■) (21) 

G {r) (X) = A - Hf-VArHjr-V- 1 

G(X) = G (r) (A)G (r - 1) (A) • • ■ G (1) (A) (22) 

then G(X) is a polynomial of A of degree r. The permut ability & implies that if 
(A a , H a ) and (Ap, Hp) are interchanged, G(X) is invariant. 
Let 

Tfij = a | 1 < a < r, l a = j} m — (m 1; • • • , m N ) (23) 

then mi + • ■ ■ + itln = r. 
Suppose 

G(A) = A r - dA^ 1 + • • ■ + (-l) r G r (24) 

denote 

C/H = itJj (Gl ) Biv ] (25) 

For the problem @ with [/ = F = 0, we have Q a = zAq,(x + Jy) + if2(\ a )t, 
where s can be t or other parameters. The matrix G(X) given by (|22] ) is a Darboux 

matrix, that is, for any solution $ of (|5|), G$ satisfies (|5|) with certain U, F, Wj, 

Xj, Zj replacing U, F, Wj, Xj, Zj. is actually the derived solution of (||). 

Owing to the choice of H a in (|17D , Z7' m l is globally defined. Here we first consider 
the generalized Q a (s) in (|T8| ) so that the localization, asymptotic properties etc. 
can be treated uniformly. 

Proposition 1 (1) If there is at most one a (1 < a < r) such that p a = 0, then 
lmw^H = o. 

(2) If p a = 0, pp = (a ^ /J), and p 1 ^ for all 7 7^ a,/?, tnen: 
(fy wnen i Q = lp, lim^oo = 0; 
(ii) when l a 7^ //3, 

limf/ir ] = for(a,b)^(l a ,lp) (26) 



and as s — > 00, 



S a/3 exp(z Im(7r a - np) + i{4>p - 4> a )s) 



W/3 -4 Q /3 cosh( Re(7r a + 7173) - o"^) + cosh( Re(7r Q - 7173) - 5^p) 

where A a p, 8 a l l, 5 a 2 }, are real constants, A a p > 0, and B a p is a complex constant. 
Moreover, if K — 1, then B a p ^ if and only if K a 7^ and up 7^ 0. 
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Proof. First suppose p a 7^ 0. By (|H|) and (|15|), 



where 



C+oo 

Q — OO 



lim H a A a H- l = Si 

p Q s-»±oo 



Aa^Af + (A Q — ^a)El a l a 



±00 



^qIk + (A a — A c 



(28 



(29) 



Ejk is an N x N matrix whose (j, k)th entry is 1 and the rest entries are zero. 
For (3 7^ a, 



(A/3 - S^°°)Hp - 



exp(g /3 (s)) 



exp(-g^( S )*)C /3 



±* 



(30) 



where 



c£ = (o,...,o,«J,o,...,o) 



A/3 — A c 
A/3 — Ac, 
A/3 -A 



K ; 3 



A/3 — A, 



A Q — A a /t/3 

A/3 — A a K*K a 



(31) 



K/3 



A Q — Aq, K*ft/} 



if lp 7^ l a 



if = L 



A/3 - A a K* a K c 

Therefore, if p a 7^ 0, the action of the limit Darboux matrix A — S^°° on Hp 
{(3 7^ a) does not change the form of Hp , but only changes the constant vector 

If K = 1, then KpK^ 7^ implies k^*R^ 0. When K > 1, this does not hold 
in general. 

Now suppose p a = 0. Without loss of generality, suppose l a = 1. Then 



I exp(7r Q ) 



1 





-exp(-7r Q X\ 




Ik 



(32) 
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By @, 

( X a A + (A a - A a ) exp(7r Q , + 7f 



A, 



(A a - A a ) exp(7r a )< \ 




V (A a - A a ) exp(7T a )K c 



A Q 
■ • • \ a AI K + (A Q - Aa)/«a< / 

(33) 

where A = exp(iT a + ix a ) + K* a K a . 

Part (1) of the proposition is derived as follows. Owing to the permutability 
of Darboux transformations, we can suppose p\ ^ 0, • • •, p r _i ^ 0, p r = 0. Then, 
as s — ► oo, tends to a diagonal matrix for a < r — 1. Considering (0), the 
limit of (G^A))^ is also diagonal, hence 



C/H = i[J, (G 1 ) Bn \ -+ 



(34) 



Now we turn to prove part (2) of the proposition. First, suppose r = 2. We 
use an alternate way of constructing Darboux matrices.0 Let 



It _ ( exp(Qa(s)) 



H* H 







A/3 — Aa 



(35) 



then 



and 



G(\) = Yl(x-x a ) l - E 

a=l V a,/3=l 



H a {r l ) a p hi 



X-Xr 



u [ri 



r 



0,(9=1 



(36) 



(37) 



Case (i): Z a = ^/3- 



7 



/I 



H n = 



exp(7r a ) 







V 
/I 



K a 



1 





(38) 



H, 



exp(7r /3 ) 



V 











t/3 







where lp — l a , then r 
Therefore, J 1-1 - 



^ where /}fc's are N x N diagonal matrices, 
where Ejk's are also N x N diagonal matrices. 



Ai A2 
. Ai A2 
^11 ^12 
.^21 ^22 
This implies that U [m] = 0. 

Case (ii): l a 7^ Suppose H a , Hp are given by (|38|) with Zg 7^ Z a . Denote 



/3 K /3 



sw? = 1 - 



4 Im X a ImXp 2 



-I* 



a/3 1 



> 



(39) 
(40) 



then, by direct computation, we have 

2i(J^ - J/J ImA Q Im\p6 a p exp(i Im(7r a - 7173)) 



A/3 — A a 



D = ^g a p cosh( Re(7r a + 7Tg) — ^) + cosh( Re(7r a — itp) — 5 2 ) 
5t = i In g aP + ^ ln(/«*ft a ^K^) + 2 In 



A/3 — A Q 



52 = Il n <^ 

2 kJK/j 



(41) 



8 



and Uau — ^ if (p, v) ^ (l a , 

When r > 2, we still use the permutability of Darboux transformations and 
suppose pi 7^ 0, • • •, p r -2 7^ 0, p r _i = p r = 0. Then, after the action of 

G( r - 2 )(A) •••£«( A), the derived F r ( ^ 2) , #( r " 2 ) has the same asymptotic form 
as H r _i, H r , provided that the constant vectors /c r _ij K r are changed to 
k^~ 2 \ Therefore, as in the case r = 2, the limit of Ui r _ lt i r has the desired form, 
and the limits of the other components of U are all zero. The proposition is proved. 

4 Localization of the solutions 

Now we consider (||), For this system, 

Q a = iX a (x + Jy) +iuo(X a )t (42) 

We consider the limit of the solution as (x, y) — > oo along a straight line x = £,+v x s, 
y = 7] + v y s (v 1 + v 2 > 0), then 

Q a = iX a (£ + Jrj) + iu(\ a )t + iX a (v x + Jv y )s (43) 

Now 

p a = Re (iX a (v x + Ji a v y )) = - Im X a (v x + Ji a v y ) (44) 

If there is at most one p a = 0, then Proposition [I] implies that £7' m ' -> as 
s — > oo. If p a = 0, Pp- = (a 7^ /3), then / a = lp since J; Q 7^ Ji if / a 7^ Z^. Hence, 

Proposition [I] also implies — > as s — ► 00. Therefore, we have 
Theorem 1 f/'" 1 ' — > as (x, y) —> oa in any directions. 

5 Asymptotic behavior of the solutions as t — > 00 

Now we use a frame (£,77) which moves in a fixed velocity (v x ,v y ), that is, let 
x = £ + t^i, V = V + v yt) then 

Q« = «A Q (£ + J77) + {iX a (v x + Jv y ) + za;(A a ))t (45) 

p a = - IhiAq,^ + J la v y ) - Im(a;j a (A a )) (46) 
Suppose that for mutually different a, /3, 7, 



det 



/ 1 a 

1 •//, ^ I / f 171 

V 1 J h ff 7 . 



where 

a a = lm(ui a (X a ))/ Im(A Q ) (4* 
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Then there are at most two p a — (a = 1, • • • , r). By Proposition [I], Ul a L ~h 
only if p a = 0, pp = 0. This leads to 

v x = J ^-y* Vy = fa^L (49) 



For UjJ? -/-> 0, a, /3 can take m 3 - and values respectively, hence there are at 

r [m] 
jk 



most m.-mjfc velocities (f x ,Vj,) such that 0. Therefore, we have 



Theorem 2 Suppose fflJQ is satisfied. Then as t — > oo ; £/ie asymptotic solution of 
Uj^ has at most m^m^ peaks whose velocities are given by (fQ) (1 a = j, lp = k). 
If a peak has velocity (v x ,v y ), then, in the coordinate £ = x — v x t, rj = y — v y t, 
\\m t -*oo U a b = for all (a, b) 7^ (j, k), and as t — > 00 

TT [m] B aj3 exp(i Re(A a - A^f + i(X a Jj - A^ J fc )r/ + z(0 a - 0^) 



jk D 

T 1 \ . 7. \„ 1 ., 

J a/3, 



D = A aP cosh( Im(A a + A^f + lm(X a Jj + \pJ k ) V + 5 { ^h < fj0 5 

+ cosh( Im(A a - A/j)£ + hn(X a Jj - A/3 Jfc)?, + %?) 

where A a p, 8^1, 8^1 are real constants, A a p > 0, and B aj3 is a complex constant, 
7 = ReA 7 (w :c + J^Vy) + Re(c^ 7 (A 7 )) (7 = a, (3) (51) 

Remark: The condition (071) implies that the velocities of the solitons are all 
different. This is true for the DSI equation. However, for the N-wave equation, all 
the solitons move in the same velocity. We will discuss this problem in the next 
section. 

Example: DSI equation 

Let n = 2, N = 2, 

Hi -0 u <-u 0) » = -™" < 52 > 

then we have 



2iJF xx + 2iUF x + i 



\u\ 2 + 2q x u x + u y \ (53) 



-u x + u y -\u\ —2q 2 
iu t = u xx + u yy + 2\u\ 2 u + 2(<7i + g 2 )w 



- = <?2,x + 92,1, = ~{\ u 



2\ 



(54) 



(FF*) D = U 9 o g ° 2 )' [^,^1 = ^ (55) 
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( |54|) is the DSI equation. 

If we construct the solution as above, then Theorem p] implies that 

{Jim] as (x, y) — > oo in any directions. If Re A a 7^ Re A/3 for a 7^ /3 and 
/ a = //3, then, theorem ^ implies that as t — > 00, w has at most m^mi peaks 
(mi + m 2 = r). From fl4q) , <r a = — 4J; Q Re A a , hence (|49|) implies that each peak 
has the velocity v x = 2 Re(A a — Ap), v v = 2 Re(A a + Ap) (l a = 1, lp = 2). This 
is the (mi, 7712) solitons.® When K — 1, these peaks do not vanish if and only if 
« a 's are all nonzero. 

Fig. 1-3 shows the solitons u^' 3 \ m' 2,3 ^ and ^t 3 ' 3 ! respectively. The parameters 
are# = 1, t = 2, Ai = 1 — 2i, A 2 = -3-i, A 3 = 2 + i, A 4 = -l+3i, A 5 = 2 + 1.5i, 
A 6 = -0.5 - 1.5i, Ci = (1, 0), C 2 = (0, 1), C 3 = (0, 2), C 4 = (0, -2), C 5 = (2, 0), 
C*6 = (—2, 0). Some techniques on drawing such figures were explained in |y|. 




Fig. 1. mI 1,3 ' of the DSI equation 




Fig. 2. J 2 ' 3 ! of the DSI equation 
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Fig. 3. J 3 ' 3 ' of the DSI equation 



6 Asymptotic solutions as the phases differences tend to 
infinity 

For the equations whose solitons move in the same speed, like the N-wave 
equation, the peaks do not separate as t — > oo. However, we can still see some 
peaks in the figures. 691 In pQ |, we obtained the asymptotic behavior of the r 2 



solitons. Here we will get the corresponding asymptotic properties of more general 
solitons obtained in this paper. We have 

Theorem 3 Let p a (a = 1. • • ■ , r) be constant real numbers satisfying 



det 



/ 1 Jl a Pa/ Im A c 

1 J h p p J Im A/? | ^ (50) 
\ 1 J; 7 pj/ Im A 7 



for mutually different a, (3, 7. Let d a be complex constant K x 1 vectors, 
K a = d a exp(p a r) and construct the Darboux transformations as above. Let 
x = £ + v x t, y = 77 + v y T, then, for any j , k with l<j,k<N,jj^ k, 

linv^oo Uj^ 7^ only if (v x ,v y ) takes specific rrijm k values. 
Proof. As in §J|, here 

Qa = i^a(t, + Jv) + iu(\ a )t + i\ a (v x + Jv y )r, (57) 

hence 

it=( eXp( ^ (r)) ) (58) 

where 

D a = (0,---,0,d aj 0,--.,0) (59) 
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Q a (r) = iX a (£ + Jrj) + iu{X a )t + (iX a (v x + Jv y ) - p a )r, (60) 
hence the real part of the coefficient of r in Q a (r) is 

p a = - lm\ a (v x + Jv y ) -p a (61) 



Condition flog) implies that there are at most two p^'s such that p a = 0. According 
to Proposition |, as r -> oo, -/->■ only if there exist p a = 0, pp = 0, a ^ (3, 
l a — j, lp — k. Therefore, the theorem is verified. 

When the condition ( fi~7| ) holds, this theorem is useless, because the evolution 
will always separate the peaks, However, when (|47|) does not hold, especially when 
it is never satisfied, this theorem reveals a fact of the separation of the peaks. 

Example: N-wave equation 

Let n = 1, u = LX where L = diag(Li, • • • , Ljv) is a constant real diagonal 
matrix such that Lj ^ Lk for j ^ k. Then, the integrability conditions (H) ~ © 
imply 

F y = JF X + UF F t = LF X + VF (62) 
[J, V] = [L, U] U t -V y + [U.V] + JV X - LU X = (63) 
U X = [J,FF*] (64) 

(|63|) is just the N-wave equation. 

Suppose is constructed as above, then Theorem |I| implies that — > 
as (x, y) — > oo in any directions. Theorem cannot be applied here. The reason 
is: the condition (|47D holds only if l a ^ lp for a ^ (3. Hence for any j, rrij = 
or 1. This implies that fliTD does not hold generally unless rrij = or 1 for all 
1 < j < TV. Therefore, we apply Theorem |3 to the previous problem. Theorem R] 
implies that if we choose {p a } such that ( jB|) is satisfied, then, for each (j,k), 

lim r ^ 00 Uj^ has at most mjm^ peaks. When K — 1, these peaks do not vanish if 
and only if k&s are all nonzero. 

Remark: Here r — >• oo means that the phase differences of different peaks tend 
to infinity. Therefore, the peaks are separated by enlarging the phase differences. 

Here are the figures describing the solutions U^- ' 1 ' 2 ^ and U^- 1 ' 1 ' 2 ^ of the 3-wave 
equation. The vertical axis is (|wi2| 2 + |wi3| 2 + |«23| 2 ) 1 ^ 4 so that all the components 
are shown in one figure. The parameters are 



J= L 



K = 1, t = 10, Ai = 1 - 2i, X 2 = -3 - i, X 3 = 2 + i, A 4 = -1 + 3z, C x = (0, 1, 0), 
C 2 = (0,0,1), C 3 = (0,0,4096), C 4 = (1,0,0). Note that for U^ X2 \ only U 23 has 
two peaks, and for Z7' 1 ' 1 ' 2 ], U%2, U13, U 23 have one, two, two peaks respectively. 
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Fig. 4. f/I 0,1 ' 2 ' of the 3-wave equation 




Fig. 5. J/! 1 ' 1 ' 2 ! of the 3-wave equation 
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